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Abstract 

Whether the a — meson (/o(600)) exists as a real particle is a long- 
standing problem in both particle physics and nuclear physics. In this 
work, we analyze the deuteron binding energy in the linear a model and 
by fitting the data, we are able to determine the range of rria and also 
investigate applicability of the linear a model for the interaction between 
hadrons in the energy region of MeV's. Our result shows that the best 
fit to the data of the deuteron binding energy and other experimental 
data about deuteron advocates a narrow range for the o"— meson mass 
as 520 < nia < 580 MeV and the concrete values depend on the input 
parameters such as the couplings. Inversely fitting the experimental 
data, our results set constraints on the couplings. The other relevant 
phenomenological parameters in the model are simultaneously obtained. 



I. Introduction 

Deuteron is the simplest bound state of nucleons, so that a study on it may provide us with 
information of the nuclear force and the fundamental principle of the strong interaction. 

Definitely nucleons are bound together by the strong interaction whose fundamental theory 
is QCD. Unfortunately, so far, we are only confident with the calculation of perturbative QCD 
at large momentum transfer i.e. at high energy scales where is sufficiently small due to 
the asymptotic freedom. By contrast, for the lower energy processes, when one evaluates the 
effects of hadronization or binding energies, the non-perturbative QCD effects dominate and the 
traditional perturbative treatments of the field theory are no longer applicable. It turns out 
that one needs to invoke some phenomenological models in the case. As a matter of fact, most 
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reasonable models express certain aspects of the real physical world, but not complete. There is 
a limited range for each model, beyond the limit, application of the model is not legitimate. On 
other side, for a certain range, all possible models arc equivalent for describing the concerned 
physics. Therefore, in some sense, various models arc parallel. This is understood as we use 
the Cornell potential, the Richardson potential or even the logarithmic potential to describe 
the bound states of heavy quarks such as J/V', T etc. and obtain close results with specific 
parameters. Of course, it by no means manifests that they are real physics, but just serve as an 
effective model and are applicable to the some concerned processes. 

For the quark-bound states, i.e. hadrons, to deal with the non-perturbative QCD effects, we 
can use, for example, the QCD sum rules, potential models, the bag model, as well as the lattice 
calculations to evaluate the spectra and other properties. It is believed that even though ag is 
large the quark-gluon picture is still valid. By contrast, for the nucleon-nucleon interaction (or 
hadron-hadron interactions when both of them are in color singlet), even though it is in principle 
due to QCD, there is no single-gluon exchange between the nucleons because they reside in 
color-singlet. Instead, just as the familiar Van der Waals force between moleculas in classical 
physics, which is induced by the electric and magnetic multipoles, the strong interaction between 
nucleons is the interaction between the chromo-electric and magnetic multipoles of nucleons. So 
far, there is no a successful way to derive an effective form of such strong interaction at the 
hadron level from the fundamental QCD theory yet, but on other side, it is also believed that 
the chiral Lagrangian is consistent with the general principles of QCD and could serve as an 
effective theory of strong interaction at hadron level [1]. An alternative version of the chiral 
Lagrangian is the linear sigma model where the cr— particle stands as an independent resonance. 
Recently, many theorists and experimentalists are searching for a which may play an important 
role in nuclear physics, because the a meson can provide reasonable middle-range nuclear force. 
However, the present data [2, 3] allow a rather wide range for both the mass and lifetime of 
(T— particle. Thus the existence of a as a real physical resonance at the hadron spectra of low 
energy is a long-standing problem in both particle physics and nuclear physics. It appears as 
/o(600) in the recent data book, but there still is an acute dispute about its existence. In this 
work, we hope to gain more information about this resonance through studying the deuteron 
structure. 

Deuteron is a familiar object in nuclear physics for many years, it may be viewed as an ideal 
lab for probing any possible mechanism of nuclear force. Many authors have recently studied 
the deuteron structure from different angles [4], and also looked for new physics, such as the 
three-nucleon force via deuteron breaking up [5] . As indicated above, we trust that some models 
are equivalent in principle for the energy range of a few MeVs, i.e. we can apply any of them. 
By fitting data, a set of concerned parameters should be obtained. 

Entem and Machleidt [6] developed an accurate NN potential based on the chiral effective 
theory. Parallel to their work, in this paper we employ the linear sigma model as an alternative 
approach to study the NN potential. The main difference from [7] is that a is an independent 
particle here and we can neglect the two-pion exchange contributions. In realistic physics world 
(T may be a resonance of certain mass and has the quantum number of two pions in an S-wave. 
The calculation now is much simpler because we do not need to calculate loop diagrams and 
deal with the renormalization. In this model, we derive the NN-effective potential and evaluate 
the spectrum, then compare it with data. In fact, besides the deuteron binding energy, there 
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are some other parameters which can be experimentally measured, such as the mixing fraction 
of the s- and d-waves, its mean charge radius etc. Therefore, when we determine the mass of 
a meson, we need also to take into account the constraints from the relevant measurements on 
the parameters. 

In fact, similar to the treatment with the chiral Lagrangian, in this work all the coefficients 
at the effective vertices are derived by fitting data of some physical processes, where all external 
particles are supposed to be on their mass-shells. Therefore, in our derivation of the effective 
potential, we need to introduce a reasonable form factor at the effective vertices which manifests 
the effects of the inner quark-gluon structure of nucleons and partly compensate the off-shell 
effects of the exchanged mesons at t-channel. Besides we also take into account the contributions 
of the light vector mesons to the potential. 

This model should be tested by evaluating the deuteron binding energy which has been 
measured with high accuracy and we wish to apply the model along with the obtained parameters 
to the discussions on the molecular states [8] which are important for understanding the hadron 
spectra, especially at middle energy ranges. 

In this work, we first derive the effective potential between nucleons. The tensor part of the 
potential results in a mixing between the s-wave and d-waves. Namely it not only contributes to 
the diagonal elements, but also to the off-diagonal ones of the hamiltonian matrix. Diagonalizing 
the matrix, we obtain the mass eigenvalues, one of them is the deuteron mass and we can argue 
that the larger one corresponds to an unstable resonance and does not exist in the physical 
world. 

This paper is organized as follows, after the introduction, we derive the formulation and in 
Sec. Ill, we present our numerical results and determine the concerned parameters, finally the 
last section is devoted to our conclusion and discussion. 



II. Formulation 



Following the traditional method [9] we derive the effective potential between proton and 
neutron. 

As discussed above we introduce a form factor to compensate the off-shell effects of the 
exchanged mesons. At each vertex, the form factor is written as [10] 

A2 - M2 
A2_g2 ' 

where A is a phenomenological parameter and its value is near 1 GeV. It is observed that as 
^ it becomes a number and if A » M^, it turns to be unity, in this case, the distance 
is infinitely large and the vertex looks like a perfect point, so the form factor is simply 1 or a 
constant. Whereas, as q'^ oo, the form factor approaches to zero, in this situation, the distance 
becomes very small, the inner structure (quark, gluon degrees of freedom) would manifest itself 
and the whole picture of hadron interaction is no longer valid, so the form factor is zero which 
cuts off the end effects. 

To derive an effective potential, one sets qq = and writes down the elastic scattering 
amplitude in the momentum space and then carries out a Fourier transformation turning the 
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amplitude into an effective potential in the configuration space. The best way to order the 
operators in the derivation is the Weyl ordering scheme [11]. 
In the linear a model, the effective Lagrangian is 



(2) 



(a) Via single-pion exchange. 
The effective vertex is 



L 



= V2g 



v2 



where ip is the wavefunction of nucleon. The scattering amplitude is 



M = [-2g pj5nnj5p + g ^75^^7571] • 



1 



mt 



(3) 



(4) 



it is noted that following the traditional treatment, the wavefunctions of nucleons in the ampli- 
tude are that of free nucleons. Setting go = and including the aforementioned form factor at 
the vertices, the potential in the momentum space has a form 



4p^ - 4(o-i • pi)(o-2 • Pi) + 2(0-2 • Pi)(o-i • q) + 2(0-2 • q)(o-i • pi) 



4m^(q^ + m^) 

- (0-1 • q)(o-2 • q) - 2(0-2 • q)(o-2 • pi) - 2(0-1 • Pi)(o-i • q) 



m: 



2x2 



A' + q' 



(5) 



where m is the mass of nTiclcon. Taking the Fourier transformation, one has the potential caused 
by the single-pion-exchange in the configuration space as 



4p?/,r(0 - 4(0-1 • Pi)(o-2 • Pi)/7r(r-) - 2i(o-2 • Pi)(o-i • r)F^(r) 



- 2iF^{r){a2 ■ r)(o-i • pi) + (0-1 • V)(o-2 • V)/^(r) + 2iF^(r)(o-2 • r)(o-2 • Pi) 
+ 2i(o-i •pi)(o-i •r)F^(r) 



(6) 



This potential is not hermitian, we need to re-order the operator p and the functions of r and 
the standard treatment is the Weyl ordering procedure [11]. The resultant effective potential is 



y-K{r)weyl 



Am? 



[/7r(r)p? + Pi/7r(r)pi + V>lU{r)] - [/7r(r)(o-i • pi)(o-2 • pi 



+ (c^i • Pi)/7r(r)(o-2 • Pi) + (0-1 • pi)(o-2 • Pi)/7r(?-)] " i[{o-2 ' Pi)(cri • r)F^(r) 
+ iSr(r)(o-i • r)(o-2 • pi)] - i[(o-i • pi)(o-2 • r)F^(r) + F„(r)(o-2 • r)(o-i • pi)] 
+ i[(o-i • pi)(o-i • r)F^(r) + F^(r)(o-i • r)(o-i • pi)] + i[(o-2 • Pi)(o-2 • r)F^{r) 

+ F,(r)(o-2 • v){<T2 ■ Pi)] + [(cri • V)(o-2 • V)/,(r)]|, 



(7) 
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where 



•""^ ' Airr 47rr SttA 

FAr) = -^U{r). (8) 

And it is the form we are going to use in the later part of the work. 

(b) Via a and p and lo exchanges. 
The effective vertices are respectively 

La = 9ipil^(T, (9) 
Lp = ^^^pV^T^r^A"'^ a =1,2, 3, (10) 

and for a;— vector meson it is 

Lu, = gj.f.u^i^-f^iljoj'^. (11) 
Neglecting some technical details, one can easily derive the corresponding effective potentials 

as 

,2 r . , [fa{r)pl + Pi/.(r)pi + p?/.(r)] 



where 



4m^ 

[VVa(r-)] , ^[Pi • rFgjr) + F^(r)r • pi] (L • S)F.(r) 



-mar -Ar / 2 _ a 2\ -Ar 



47rr 47rr SttA 



and 

Vp{r) 



^{4m^ /,(r) + [V^ Ur)] - [/p(^)Pi + Pi + Pi/p(0] + 2F,(r)L • S 

5 

- z[pi • ri?;,(r) + -Fp(r)r • pi] + -[fp{r){cri ■ pi)((T2 • Pi) + {en ■ pi)/p(r)((T2 • Pi) 
+ (cTi • pi)(o-2 • Pi)/p(r)] + y [Fp(r)(fTi • r)(o-2 • Pi) + (0-2 • Pi)(cri • r)Fp(r)] 
+ ^[^pW('T2 • r)((Ti • pi) + ((Ti • pi)(o-2 • r)F,(r)] - [{a, ■ V)(o-2 • V)/,(r)] 
+ ^[(o-i X 0-2) • pi(cri • r)Fo(r) + Fp{r){ai ■ r){(Ti x 0-2) • pi] 
+ 7;[Fp{r){cri X 0-2) • r(cr2 • Pi) + {(T2 ■ Pi)((ri x ■ rFp{r)] 



- 2^Ppir)i(^i X 0-2) • r(<7i • pi) + (o-i • pi)((Ti x (T2) • rFp{r)] 

- X 0-2) • Pi(<T2 • r)Fp(r) + Fp{r){cr2 ■ r)((Ti x (T2) • Pi]|, 



(14) 



where 



fp{r) 



-Ar 



Airr 

id_ 

r dr 



Airr 



+ 



SttA 



-Ar 



Z-Jpir). 



(15) 



The potential induced by exchanging u vector meson is in close analogy to Vp, with only nip 
being replaced by rriuj. We neglect the contributions by exchanging heavier mesons. 

The total effective potential is the sum of these potentials induced respectively via exchanging 
TT, (7 and p, uj mesons. 

Veff{r) = 14 (r) + K(r) + Vp{r) + V^{r). 

After simple manipulations, we can write the Veffi^) as a sum containing a central potential 
term Vc{r), a spin-dependent central potential part Vs{r) <T\ ■ <T2, a spin-orbit coupling part 
yLs{f) L • S, and a tensor interaction part VT{r) Si2- The tensor operator 5*12 = [3(<ti ■ f)(<T2 • 
f ) — <Ti • 0-2] is known to admix the s-wave and d-wave states, so we separate out the tensor part 
of the potential as 



VT{r)S: 



12 



-Ar^3 2 



+ 



+ 



64m^7r 

-Jp-mpr 2 



+ 



''Amlg 



P^nnP 



nnP 



64m^7r 



32m2 



+ 



NNP 



nnP 



32m2 



7rr^ 



3e-^'-A2/ 



dAm^TTr 



-Ar 2 2 



'UJiJ l^l^W 



dAm^TTr 



P^NNP _j_ ^ '"'pi' NNP 



3e 



64m^7r 

■^^A^L,,,, 3e 



32m27rr2 



+ 



64m^7r 

2 



32m2 
Se-^'-A^^^ 



o -ma,r 2 

32m27rr^ 

,2 



32m27rr 



32m27rr2 
'S'12, 



7rr 



Trr 



(16) 



and treat it as a perturbation later. 

Now, the Schrodinger equations for s-wave state (^Si) and d-wave state (^Di) have the 
following forms 



2/x 



+ Vc{r) + Vs{r) 



*,(r) = S,*,(r), 



+ Vc{r) + Vsir) - 3VLs{r) - 2VT(r) 



*d(r) = £;d*d(r), 



(17) 



(18) 



where /x is the reduced mass, here it is ^Mp = ^M„. Solving these equations, we obtain the eigen- 
energies and eigen-wavefunctions of the s-wave and d-wave states. Then the tensor part, which 
results in a mixing fraction between the s- and d-waves, would be taken in as a perturbation. 
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We have a modified hamiltonian as 

H,^,= [f. ^J. (19) 

Here Eg and are the eigenvalues of s- and d-waves obtained by solving the Schrodinger 
equations (17) and (18), A is defined as 

A=< s\VTir)Si2\d > . 

Diagonalizing HjuU. we obtain the real eigenvalues and identify the lower one to be the binding 
energy of deuteron (we will discuss this issue in the last section again) , at the same time we can 
determine the fraction of the d-wave in the eigenfunction 

I* >= ci\s > +C2\d > . 

III. The numerical results 

The input parameters include the masses = 0.139 GeV, nip = 0.77 GeV, rriuj = 0.783 
GeV, nip « nin = 0.938 GeV, and the couplings g^^j^a = 9nnt^ — 13.5, ^jvatp — 3-25 [12]. In the 
exact SU(3) limit, by a simple quark counting, g^^j^oj = "^gNNp-: but our later numerical results 
show that this relation does not lead to a satisfactory solution. For a comparison, we will take 
three difi'erent gf^^^uj values which are gj^j^p-, ^g^NP Sgj^a^p respectively and present the results 
in tables 1,2 and 3. 

Even though we do not have direct information about gNNa, gNNa = gNNn is imposed by 
the linear sigma model. A is a free parameter in the form factor, it is generally believed that 
its value is in a range around 1 GcV. Therefore, one can vary it within a reasonable region 
around 1 GeV. Here we take it as a free parameter close to 1 GeV and adjust it to fit the data of 
deuteron. When we carry out the numerical computations, by imposing the experimental data 
as constraints which include the binding energy of deuteron, a mixing between s- and d-waves 
and the charge radius etc., it is found that a unique A— value for both s- and d-waves cannot 
lead to a reasonable solution set which can meet the data. Therefore we take separate values 
for s- and d-waves, more discussions will be given in the last section. 

The main task of this work is to find the mass of cr— particle, it stands as free a parameter 
in the present calculation and we will scan the possible region and finally determine its value 
by fitting data. As well known, a— meson is a rather wide resonance, it is so far not well 
experimentally determined yet, the data allow a large range as the mass being 400 ~ 1200 MeV 
and width being 600 ~ 1000 MeV[2]. For deriving the effective potential the width is irrelevant, 
because the t-channel is a deep space-like region, so the width is not involved, thus in this work 
we cannot determine the width from the static properties of deuteron. If one wishes to study 
the width, he must calculate the production rate or decay (lifetime) of deuteron which would 
be more difficult and we will deal with it in our later work. 

The experimental value of the binding energy of deuteron is [15] 

^exp ^ _2.224575 MeV. 

Our strategy is to fit the binding energy which is very accurately measured and serves as a 
rigorous criterion for our computation. Then we achieve a reasonable mixing fraction between 
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the s- and d-waves which has also been experimentally confirmed. As aforementioned, to fulfil 
all the criteria, a unique A— value is not enough, thus we set A^ and A,^ for s- and d-waves 
respectively when we solve the Schrodinger equations. We present our numerical results in 
Table 1. 



ma 

(xlO-^GeV) 


A. 
(GcV) 


Ad 
(GeV) 


Es 
(MeV) 


Ed 
(MeV) 


A 

(MeV) 


\Cs> 


r 
(fm) 


4.70 


0.59 


1.35 


-1.97 


-0.29 


0.71 


0.13 


3.54 


4.80 


0.60 


1.40 


-2.00 


-0.36 


0.64 


0.12 


3.47 


4.90 


O.Gl 


1.16 


-2.01 


-O.-IO 


0.58 


0.10 


3.39 


5.00 


0.62 


1.52 


-2.08 


-0.43 


0.52 


0.08 


3.31 


5.10 


0.63 


1.58 


-2.11 


-0.44 


0.46 


0.07 


3.23 


5.20 


0.64 


1.65 


-2.13 


-0.48 


0.40 


0.05 


3.16 


5.30 


0.65 


1.71 


-2.15 


-0.51 


0.35 


0.04 


3.08 


5.40 


0.66 


1.78 


-2.17 


-0.54 


0.29 


0.03 


3.01 


5.50 


0.67 


1.84 


-2.19 


-0.55 


0.25 


0.02 


2.94 


5.60 


0.69 


1.91 


-2.22 


-0.58 


0.20 


0.01 


2.88 



Table 1. gNNir = QNNa = 13.5,57ViVp = 3.25,5jvAra; = 9NNp All the other parameters are 
obtained by fitting the binding energy E^^^ = —2.22 GeV 

With another set of couplings, we have 



(xlO-^GeV) 


A. 
(GeV) 


Ad 
(GeV) 


Es 
(MeV) 


Ed 
(MeV) 


A 

(MeV) 


(8)= 


r 

(fm) 


5.10 


0.64 


1.75 


-1.70 


-0.30 


1.00 


0.27 


3.54 


5.20 


0.65 


1.84 


-1.80 


-0.33 


0.90 


0.22 


3.45 


5.30 


0.66 


1.92 


-1.89 


-0.36 


0.79 


0.18 


3.36 


5.40 


0.67 


2.02 


-1.97 


-0.40 


0.68 


0.14 


3.28 


5.50 


0.68 


2.11 


-2.04 


-0.43 


0.58 


0.10 


3.19 


5.60 


0.69 


2.21 


-2.09 


-0.46 


0.48 


0.08 


3.10 


5.70 


0.70 


2.30 


-2.13 


-0.50 


0.39 


0.05 


3.02 


5.80 


0.71 


2.40 


-2.17 


-0.51 


0.31 


0.03 


2.92 



Table 2. gNNn = QNNa = 13.5,(7Arjvp = 3.25, gNNuj = ^g^Np AH the other parameters are 
obtained by fitting the binding energy E^^^ = —2.22 GeV 

Then for a comparison we give another set of results which is listed in table 3. 



8 



(xlO-^GeV) 


As 
(GeV) 


Ad 
(GeV) 


Es 
(MeV) 


Ed 
(MeV) 


A 
(MeV) 




r 
(fm) 


3.50 


0.53 


0.79 


-2.16 


-0.44 


3.81 


0.64 


4.45 


4.00 


0.57 


1.03 


-2.08 


-0.40 


3.46 


0.62 


3.97 


4.50 


0.61 


1.46 


-1.92 


-0.25 


2.72 


0.55 


3.67 


5.00 


0.65 


2.14 


-1.84 


-0.19 


1.86 


0.43 


3.44 


5.50 


0.69 


3.07 


-1.68 


-0.09 


1.02 


0.24 


3.36 



Table 3. gNNn = gNNa = 13.5,(7ArArp = 3.25,gis[Nu) = '^dNNp All the other parameters are 
obtained by fitting the binding energy E^^^ = —2.22 GeV. 

Obviously, even though all the three sets of couplings can lead to results which coincide with 
the data on the binding energy, when we consider other experimental constraints, especially the 
s-wave and d-wave mixing fraction of about 4%, the set given in table 3 is not favored. 

It is noted that to fit the binding energy of deuteron and consider the constraint of the s- and 
d-wave mixing, the estimated mass of a meson deviates for different parameter sets by about 50 
MeV. We will further discuss this issue in the last section. And we have also noted that for each 
parameter set, beyond the range of rrifj^ we cannot find any reasonable solution for the binding 
energy and the mixing parameter {Cd/Cgf'- 

IV. Discussions and conclusion 

Nowadays, the hadron structure is still a hot topic in particle physics, we wonder if the 
observed hadron spectra from light to heavy can be fully interpreted in the quark model. Among 
all the light hadrons, the scalar a (0^, /o(600)) may take the most noticeable position. Dispute 
about its existence as a physical resonance has been on for a long time. A recent observation of 
the narrow resonance (2.32 GeV/c^) by the BABAR collaboration [19] has triggered a new tide 
to discuss the whole family of the 0"*" particles [20] where a is the leading one. One of the most 
striking questions is whether they can be attribiitcd to the building of qq system or some of them 
may be molecular states. Indeed, the measured resonance peaks /o(980), ao(980) are interpreted 
as KK molecula [8], Rujula, Georgi and Glashow suggest that ^'(4.04) is a D*D* molecula [14] 
etc. To further investigate such molccTdar states, one needs to know the interaction between 
the meson-constituents in the molecular states and derive the effective potential. In that case, 
the a particle may be the main agent to mediate the interaction as well as pion and other light 
vector mesons. Thus to test the properties of a becomes the most important task in the whole 
study. 

There are many ways to study a, among them the structure of deuteron provides us with 
direct information about the a particle. Deuteron is the simplest bound state of nucleons, in 
general, it is a "molecular" state of two fermions (proton and neutron). Because its spectrum 
is accurately measured, it is an ideal probe for the applicable mechanism and theoretical frame- 
work. 

In this work, we employ the linear a model where the cr is a real physical particle whereas 
in the chiral Lagrangian approach its role is replaced by the two-pion-exchange picture. At the 
effective vertices, we retain a form factor with a free parameter A which is determined by fitting 
data. 
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When we consider the constraint from the experimental observation that a certain fraction 
of about 4% d-wave mixes with the s-wave, we fit the binding energy and the results are listed in 
Table 1. As aforementioned, a unique A value cannot lead to a solution for the binding energy 
and mixing of around 4%. Thus wc take two parameters A^ and corresponding to s- and 
d- waves respectively. In fact, A represents the inner structure of the constituents of deuteron 
and the effective interaction between them. For s- and d-waves, the orbital angular momenta 
are different, and the momentum distributions of the constituents in the s- and d-states behave 
differently, namely the effective binding could deviate from each other. The different A— values 
for s- and d-waves would manifest the difference in effective interaction for different angular 
momenta states. 

The effective coupling gj^j^a = 9NN-n- = 13.5 is taken in our numerical computations, but other 
possibilities are also noted that in literature, a different value of the coupling is also commonly 
used as ^jv^a = g^N-^ — 10.1[16]. To investigate the influence of different couplings, we repeat 
all the numerical computations with the coupling being replace by gj^j^a = ffjvivTr = 10.1, the 
results are not in a good shape, namely no reasonable solutions for both the binding energy and 
mixing fraction can be reached. Thus this model does not advocate smaller coupling gj^j^a = 
9nnt^ = 10.1. 

It is worth noticing that the results somewhat depend on the coupling gj^j^u- Indeed, both 
the p and u mesons provide an effective short-range expelling interaction between nucleons, by 
the SU(3) symmetry, ^j^j^,/, = and gj^^^ui ~ "^Onnp [12]) if a simple quark counting rule is used. 
But the relation does not need to be taken seriously, because the SU(3) symmetry is obviously 
broken. In fact, in the literature, quite different relations between ^^vjvu; and gj^ffp have been 
adopted. One would rather use the loose relation as S'^vatw = (^9nnp where a is treated as a free 
parameter and is to be determined by fitting data. In the text, we adopt three different a values 
to be 1, 2 and 3 respectively and present corresponding results in the tables. We observe that 
all the results corresponding to the three relations can lead to solutions which perfectly coincide 
with the data of the deuteron binding energy, but while considering the s-wave and d-wave 
mixing fraction {cd/csl"^ as a constraint, one can immediately notice that for larger a values the 
fraction becomes unbearably large. It indicates that smaller a values are more favored, and it 
may also manifest the SU(3) breaking in this case or in other words, for an effective theory for 
the NN interaction at very low energy (near zero), g'^v^oj ~ 9nnp- 

As a conclusion of this work, the data on deuteron advocates the mass range of the a particle 
around 520 to 580 MeV which is consistent with the present data and the value commonly used 
in the phenomenological nuclear physics. Meanwhile we obtain a fraction of the d-wave in the 
total wave function as about 4~ 10% and it also coincides with the observation. One question 
might be raised that if deuteron is a mixed state of the s- and d-waves with the lower eigen- 
energy, where is the bound state of the higher energy. Since deuteron is the unique stable nucleus 
consisting of a proton and a neutron, no other bound state has been experimentally observed, 
it implies that the bound state of higher eigen-energy is unstable and easy to dissociate, so does 
not exist in the nature. Indeed, our calculation shows that as the binding energy of deuteron is 
—2.22 MeV, another cigcn-cncrgy is very close to zero. 

It has been reported that the newly measured value of the mass and width of cr— meson as 

m^ = 390l|gMeV and T,, = 282^^^ MeV[17]. 
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Huo et al. recently also discussed on the a properties in J ftp decays [18]. Their value of the 
mass of cr— meson might be a bit smaller than the generally expected. Instead, in this work 
we use the deuteron data to investigate the mass. Because of the measurement errors, we take 
several mg- values and recalculate Ei,. 

In general, our result indicates that the linear a model may apply to determine phenomeno- 
logical quantities of deuteron and inversely by studying these quantities, one can achieve further 
information about this long-expected particle /o(600). As many theorists and experimentalists 
concern its existence and search for it at various experimental facilities and processes, we be- 
lieve that its existence will be confirmed soon and its properties including mass and lifetime 
will be fixed with much higher accuracy in near future. Then we will be able to further study 
our mechanism and the parameter A. These results can be compared with other data of wider 
energy ranges for getting a better insight into the strong interactions. That achievements would 
enrich our knowledge about the interaction mediated by the a meson and then we will apply 
it to other areas of high energy physics and nuclear physics. However, from the other side, we 
find that the results depend on the adopted parameters which are obtained by fitting the data 
of nucleon scattering and given in different works of literature. Further study on a meson would 
provide more definite information about its mass and then we can continue to investigate the 
linear a model. 
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